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Abstract 

We calculate the van der Waals friction between two semi-infinite 
solids in normal relative motion and find a drastic difference in com- 
parison with the parallel relative motion. The case of the good con- 
ductors is investigated in details both within the local optic approxi- 
mation, and using a non-local optic dielectric approach. We show that 
the friction may increase by many order of magnitude when the sur- 
faces are covered by adsorbates, or can support low-frequency surface 
plasmons. In this case the friction is determined by resonant photon 
tunneling between adsorbate vibrational modes, or surface plasmon 
modes. The theory is compared to atomic force microscope experi- 
mental data. 

1 Introduction 

A great deal of attention has been devoted to non-contact friction between 
nanostructures, including, for example, the frictional drag force between two- 
dimensional quantum wells PHI21IIII ; and the friction force between an atomic 
force microscope tip and a substrate |H El El HI IE] • 

In non-contact friction the bodies are separated by a potential barrier 
thick enough to prevent electrons or other particles with a finite rest mass 
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from tunneling across it, but allowing interaction via the long-range elec- 
tromagnetic field, which is always present in the gap between bodies. The 
presence of inhomogeneous tip-sample electric fields is difficult to avoid, even 
under the best experimental conditions |6 ; . For example, even if both the tip 
and the sample were metallic single crystals, the tip would still have corners 
present and more than one crystallographic plane exposed. The presence of 
atomic steps, adsorbates, and other defects will also contribute to the inho- 
mogeneous electric field. The electric field can be easily changed by applying 
a voltage between the tip and the sample. 

The electromagnetic field can also be created by the fluctuating current 
density, due to thermal and quantum fluctuations inside the solids. This 
fluctuating electromagnetic field is always present close to the surface of any 
body, and consist partly of traveling waves and partly of evanescent waves 
which decay exponentially with the distance away from the surface of the 
body. The fluctuating electromagnetic field originating from the fluctuating 
current density inside the bodies gives rise to the well-known long-range at- 
tractive van der Waals interaction between two bodies jHj. If the bodies are 
in relative motion, the same fluctuating electromagnetic field will give rise to 
a friction which is frequently named as the van der Waals friction. Van der 
Waals friction can be considered as mediated by photon exchange between 
the bodies: One body emit a photon, and the other absorbs it, thus trans- 
ferring momentum between the bodies, resulting in a friction force. At large 
distances between the bodies, the main contribution to friction comes from 
photon exchange, corresponding to the propagating electromagnetic waves. 
However this contribution is very small because the photons corresponding to 
propagating waves carry very small momentum, no larger than ksT/ch. The 
photons, corresponding to the evanescent electromagnetic waves, carry the 
momentum q ~ d^ 1 . Thus for distances d between two bodies smaller char- 
acteristic distance dx = hc/hsT, which depends on temperature (at room 
temperature dx ~ 10 5 A), the main contribution to friction comes from the 
evanescent electromagnetic field. In analogy with electron tunneling, this 
mechanism of momentum transfer can be considered as associated with the 
photon tunneling. 

Although the dissipation of energy connected with the non-contact fric- 
tion always is of electromagnetic origin, the detailed mechanism is not to- 
tally clear, since there are several different mechanisms of energy dissipation 
connected with the electromagnetic interaction between bodies. First, the 
electromagnetic field from one body will penetrate into the other body, and 
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induce an electric current. In this case friction is due to Ohmic losses in- 
side the bodies. The fluctuating electromagnetic field can also excite the 
vibrations of the adsorbates or other surface localized modes, e.g. surface 
plasmons and polaritons. In this case friction is due to energy relaxation of 
the surface modes. Another contribution to friction from the electromagnetic 
field is associated with the time- dependent stress acting on the surface of the 
bodies. This stress can excite acoustic waves, or induce time-dependent de- 
formations which may result in a temperature gradient. It can also induce 
motion of defects either in the bulk, or on the surface of the bodies. The 
contribution to friction due to non-adiabatic heat flow, or motion of defects, 
is usually denoted as internal friction. 

It is very worthwhile to get a better understanding of different mecha- 
nisms of non-contact friction because of it practical importance for ultrasen- 
sitive force detection experiments. This is because the ability to detect small 
forces is inextricably linked to friction via the fluctuation-dissipation theo- 
rem. For example, the detection of single spins by magnetic resonance force 
microscopy, which has been proposed for three-dimensional atomic imag- 
ing 14 and quantum computation [12], will require force fluctuations to 
be reduced to unprecedented levels. In addition, the search for quantum 
gravitation effects at short length scale ^B] and future measurements of the 
dynamical Casimir forces may eventually be limited by non-contact fric- 
tion effects. 

Recently Gotsmann and Fuchs jH] observed long-range non-contact fric- 
tion between an aluminum tip and a gold (111) surface. The friction force F 
acting on the tip is proportional to the velocity v , F = Tv. For motion of the 
tip normal to the surface the friction coefficient T(d) = b- d~ 3 , where d is the 
tip-sample spacing and b = (8.0±|;|) x l(r 35 Nsm 2 0. Later Stipe et.al.fl 
observed non-contact friction effect between a gold surface and a gold-coated 
cantilever as a function of tip-sample spacing d, the temperature T, and the 
bias voltage V. For vibration of the tip parallel to the surface they found 
T(d) = a(T)(V 2 + V£)/d n , where n = 1.3 ± 0.2, and V ~ 0.2 V. At 295K, 
for the spacing d = 100A they found T = 1.5 x 10 _13 kgs _1 , which is ~500 
times smaller that reported in Ref. j3] at the same distance using a parallel 
cantilever configuration. 

In a recent Letter, Dorofeev et. al. jl] claim that a the non-contact friction 
effect observed in jU |3] is due to Ohmic losses mediated by the fluctuating 
electromagnetic field . This result is controversial, however, since the van der 
Waals friction has been shown [TUllTHir2*l lT3] to be many orders of magnitude 
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smaller than the friction observed by Dorofeev et.al. Presently, the origin of 
the difference in magnitude and distance dependence of the long-range non- 
contact friction effect observed in [5^ and P is not well understood. 

In order to improve the basic understanding of non-contact friction, we 
present new results for van der Waals friction. In JU] we developed a the- 
ory of van der Waals friction for surfaces in parallel relative motion. Here 
we generalize the theory to include also the case when the surfaces are in 
normal relative motion, and we show that there is drastic difference between 
these two cases. Thus, for normal relative motion of clean good conductor 
surfaces, the friction is many orders of magnitude larger than for parallel 
relative motion, but still smaller than observed experimentally. Another 
enhancement mechanism of the non-contact friction can be connected with 
resonant photon tunneling between states localized on the different surfaces. 
Recently it was discovered that resonant photon tunneling between surface 
plasmon modes give rise to extraordinary enhancement of the optical trans- 
mission through sub- wavelength hole arrays [IS] . The same surface modes 
enhancement can be expected for van der Waals friction if the frequency 
of these modes is sufficiently low to be excited by thermal radiation. At 
room temperature only the modes with frequencies below ~ 10 13 s _1 can be 
excited. For normal metals surface plasmons have much too high frequen- 
cies; at thermal frequencies the dielectric function of normal metals becomes 
nearly purely imaginary, which exclude surface plasmon enhancement of the 
van der Waals friction for good conductors. However surface plasmons for 
semiconductors are characterized by much smaller frequencies and damping 
constants, and they can give an important contribution to van der Waals 
friction. Other surface modes which can be excited by thermal radiation are 
adsorbate vibrational modes. Especially for parallel vibrations these modes 
may have very low frequencies. 

All information about the long-range electromagnetic interaction between 
two non-contacting bodies is, in principle, contained in the reflection factors 
of the electromagnetic field. At present time very little is known about the 
reflection factors for large wave vectors and for extremely small frequencies. 
In the calculations of the reflection factors one must take into account the 
non-local response of the electron gas on the external electromagnetic field. 
There are two correlation length which determine this nonlocal response. 
First is the skin depth, which determines the long-range length scale of the 
nonlocality in the volume, and the second is the screening length. The lat- 
ter length scale determines the short range correlation length for non-local 
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response in the surface region. In our previous calculations of the Van der 
Waals friction (TU1 HU E21 EE] we mostly considered good conductors. In this 
case it was shown that the important contribution comes from the non-local 
optic effects in the surface region. However it was shown that the Van der 
Waals friction becomes much larger for high resistivity material, for which 
the volume contribution from non-local effects is also important. It is easy 
to see that within local optic approximation the Van der Waals friction di- 
verge when the conductivity of materials tend to zero. This means that the 
local optic approximation breaks down for high-resistivity materials. This 
situation is completely different from the heat transfer between bodies via 
photon tunneling [TTH IT2] , where the heat flux is maximal at conductivities 
corresponding to semi-metals. In order to clarify the situation we study the 
dependence of the van der Waals friction on the dielectric properties of the 
materials within the non-local dielectric approach, which was proposed some 
years ago for the investigation of the anomalous skin effects [20J. 

2 Calculation of the fluctuating electromag- 
netic field 

We consider two semi-infinite metals 1 and 2 having parallel flat surfaces. We 
introduce a coordinate system with xy plane in the surface of body 1 , and the 
z axis along the upward normal. The surface of body 2 is located at z — d, 
performing small amplitude vibrations along the z axes with displacement 
coordinate u z (t) = UQe~ luJot . Since the system is translation invariant in the 
x = (x, y) plane, the electromagnetic field can be represented by the Fourier 
integral 



where E and B are the electric and magnetic induction field, and q is the two- 
dimensional wave vector in (x, y) plane. After the Fourier transformation it is 
convenient to choose the coordinate axis in the (x,y) plane along the vectors 
q and n=[zxq]. In the vacuum gap between the bodies the electric field 
E(q, u,z), and the magnetic induction field B(q, u),z), can, to the linear 
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order in the vibrational coordinate, be written in the form 

E(q, u, z) = ((w Q e ipz + v e- ipz ) + (w ie ip+z + v ie - ip+z ) e - ^ ') e~ iuJt (3) 



B(q, u, z) = c 



u 

+ — - — ([k; x v l ]e~ ip+z + [k+ x w ie ip+z }) e^"'"' 

where = q ± zp, p = ((lu/c) 2 — q 2 ) 1 ^ 2 , p + = p(uo + loq), k+ = k + (cj + 
uq), k^ = k~(u + uo). At the surfaces of the bodies the amplitude of the 
outgoing electromagnetic wave must be equal to the amplitude of the reflected 
wave plus the amplitude of the radiated wave. It is convenient to decompose 
the electromagnetic field into the p- and s - polarized electromagnetic waves. 
For p-polarized electromagnetic waves the electric field is in the incident plane 
determined by the vectors q and n, and for s- polarized electromagnetic waves 
the electric field is normal to the incident Thus the boundary conditions for 
the electromagnetic field at z = can be written in the form 

Woz(y) = Rip(s)(uj)v 0z (y) + E f lz(y) (u) (5) 

Wlz(y) = Rlp(s)(u + ^o)Vi z (y) (6) 

where Ri p ^(u>) is the reflection factor for surface 1 for p(s) - polarized elec- 
tromagnetic field, and where E( z ^ (u) are the components of the fluctuating 
electric field outside the surface 1 in the absence of the body 2. The boundary 
condition at the surface of the body 2 must be written in the reference frame 
where the body 2 is at rest. The electric field in this reference frame is de- 
termined by a Lorentz transformation. Performing a Lorentz transformation 
of the electric field to linear order in uq gives 

E , = E _ iu u(t) [e z x B] ^ 

For the p-polarized electromagnetic waves the second term in (|7J) is of the 
order of magnitude ujqUqoj /pc 2 relative to the first one and can be neglected 
for the most practical cases. However, for the s-polarized electromagnetic 
waves the second term is of the order of magnitude uj u p/uj and can be of 
the same order of magnitude as the first term. In the rest frame of body 2 



- (V x v ]e- ipz + [k+ x w ]e ip ' 

e" iwt (4) 
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there is also mixture of s- and p— polarized electromagnetic waves. In ^U] it 
was was shown that this gives contribution of the order {ujqUq/c) 2 and thus 
can be neglected. After performing Lorentz transformation to linear order in 
ujq and u we get Vq = v , Wq = w 

v 'iz(x) = v iz(x) - ipu v 0z ( x y, w' lz{x) = wi z ( x) + ipu w 0z ( x y, 
. UJ + tu . . uj + lu . 

W ly = W ly ^ ~ ipU W y] % = V ly ipU V 0y 

The boundary conditions for the electromagnetic field at z = d + u(t) in the 
rest frame of body 2 can be written in the form 

v 0z{ y) = e 2ipd R Ms) (u)w 0z{y) + e^ d E{ z{y) (uj) (8) 
vu - ipu v 0z = e 2ip+d R£ p (w lz + ipu w 0z ) (9) 

(u + UJ )v 0y 2ipdo+ ( , ■ (uJ + Uo)w y\ . . 

v ly - ipu = e p R^ s \ Wiy + tpu y - (10) 

uj V uj J 

where R2 P ( s )(u) is the reflection factor for surface 2 for p(s) - polarized elec- 
tromagnetic field, and where E^ z ^ (uj) are the components of the fluctuating 
electric field outside the surface 1 in the absence of the body 1. From ()5I6|) 
and ABE! we get 

RMs)EL {y) e^ d + E[ z[y) 

W 0z(y) = ^ (UJ 

e 2tpd RM s) E{ z{y) + E{ z{y) e^ d 

V 0z (y) = ^ (12) 

. (e 2 ^ d R f 2p + e 2i P +d Rt p )E{ z + (1 + e 2 ^ d Rt v R x )E{ z e^ d 
v lz = ipu v - — — — (13) 

. uj + uj (e 2ipd R 2s + e 2i P +d R+)E[ y + (1 + e*r +d R+ R^E^ 

Vly = IPUQ ■ — 



uj A S A+ 

(14) 

Wlz(y) = Rt p[s) Vlz(y) (15) 

where R^ s = R p{s) (uj+uj ), A p{s) = l-e 2i P d R 2p ( s )Rip( s ), and A+ (s) = A p(s) (uj+ 
uj ). Other components of the fluctuating electromagnetic field can be found 
from the transversality conditions 

qw x + pw z = 0, qv x - pv z = (16) 
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The fundamental characteristic of the fluctuating electromagnetic field is 
the correlation function, determining the average product of components 
E*(q, u). Accordingly to the general theory of the fluctuating electromag- 
netic field (see for a example J2|) these correlation function are given by 

< |^(q,^)| 2 >= (n(u>) + £) [(p+p*)(1-\Rs\ 2 ) + ( P -P*)(R:-Rs)} 

(17) 

< |£/(q,u;)| 2 >= |J (n(fj) + i) [(p+P*)(l - \R V \ 2 ) + (p- P *)(R; - R,)) 

P (18) 
where < ... > denote statistical average over the random field, and where the 
Bose-Einstein factor 



We note that p is real for q < u/c (propagating waves), and purely 
imaginary for q > u/c (evanescent waves). Thus for q < u/c and q > u/c 
the correlation functions are determined by the first and the second terms in 
Eqs. (JT7j)and (JTHJ), respectively. 



3 Calculation of the friction force between 
two semi-infinite bodies in normal relative 
motion 

The frictional stress a which act on the surfaces of the two bodies can be 
obtained from zz— component of the Maxwell stress tensor , evaluated 
at z = 0: 

° zz = An Jo duJ J (27r)2[ < \ Ez ^ u ' z ^ 2 > + < l^(q>^)| 2 > 

- < (^(q,^,^)! 2 > - < \E y (q,u,z)\ 2 > 

- < \B x (ci,u;,z)\ 2 >-< \B y (q,uj,z)\ 2 > (19) 

To linear order in the vibrational coordinate u(t) and the frequency ujq , the 
stress acting on the surface 1 can be written in the form 

d 

cr Z z = cr 0zz (d) + u(t)—a 0zz (d) + iuoj^t) (20) 
8 



Here the first term determines the conservative van der Waals stress and the 
second term is the adiabatic change of the conservative van der Waals stress 
during vibration. The last term determines the frictional stress with friction 
coefficient 7 ± . For normal relative motion (see Appendix A) we obtain the 
friction coefficient / y ± = Y± d + 7^ mn , where the contribution to the friction 
coefficient from the propagating electromagnetic waves is given by 

uj/c 

dq qp 2 

x[(l - |iVR 2p | 2 ) 2 + |(1 - \R lp \ 2 )R 2p e ipd 




+ (1 - W ^ r1,| M^ + ^ 4 (21) 

and where the contribution to the friction from the evanescent electromag- 
netic waves is given by 

7T Z Jw/c \ OUJ J Juj/c 

x [(JmR lp + e~ 2kd \R lp \ 2 lmR 2p )(lmR 2p + e~ 2kd \R 2p \ 2 lmR lp ) 

+ e -^(I m(Wfl — l— l + lp ^ s] , (22) 

where k = \p\. The symbol [p — > s] in Eqs. (121 J) and (|2*2*jl denotes the 
term which is obtained from the first one by replacement of the reflection 
factors Rp(uj), for p— polarized waves, by the reflection factors R s (w) for 
s— polarized waves. The friction coefficient for two flat surfaces in paral- 
lel relative motion was obtained by us before, JU] and can be written as 
7|l = 7|j ad + 7jj mn , where the contribution to the friction coefficient from the 
propagating electromagnetic waves is given by 

P, pOO I fir, \ rUl/c 



|1 - e 2t P d R lp R 2p \ 
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and where the contribution to the friction from the evanescent electromag- 
netic waves is given by 



rad 




(-¥) r ^ 

\ OOJ ) Ju/c 



-2kd 



1 



2 + \P 4 



(24) 



x ImRiplmR 



'2p 



1 - e^ kd R lp R 2p 



There is a principal difference between the friction coefficient for normal and 
parallel relative motion, related to the denominator in the formulas for the 
friction coefficient. The resonant condition corresponds to the case when the 
denominator of the integrand in Eqs. (j21H24|) . which is due to multiple scat- 
tering of the evanescent electromagnetic waves from the opposite surfaces, 
is small. For two identical surfaces and Ri « 1 < R r , where Ri and R r 
are the imaginary and real part, respectively, this corresponds to the reso- 
nant condition i?^exp(— 2kd) ~ 1. At resonance the integrand in Eqs. (|2*Hj) 
and (J2"4"j) has a large factor ~ ^/R 2 , in sharp contrast to the case of parallel 
relative motion, where there is no such enhancement factor. The resonance 
condition can be fulfilled even for the case when exp(— 2kd) << 1 because for 
evanescent electromagnetic waves there is no restriction on the magnitude of 
the real part or the modulus of R. This open up the possibility of resonant 
denominators for R 2 r » 1. 

To estimate the friction coefficient T for an atomic force microscope tip 
we can use an approximate formula [2H 122] 



where it is assumed that the tip has cylinder symmetry. Here z(p) denotes the 
tip - surface distance as a function of the distance p from the tip symmetry 
axis, and the friction coefficient ^(z(p) is determined by the expressions for 
the flat surfaces. This scheme was proposed in [21J for the calculation of 
the conservative van der Waals interaction. The error of these scheme is not 
larger than 5-10% in practice in an atomic force microscopy experiment, and 
25% in a worst case satiation [22J. Although this scheme was proposed for 
the conservative van der Waals interaction, we assume that the same scheme 
is also valid for the calculation of the van der Waals friction. We assume that 
the tip has a paraboloid shape given [in cylindrical coordinates (z, p)} by the 
formula:^ = d + p 2 /2R, where d is the distance between the tip and the flat 




(25) 
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surface, and where R is the radius of curvature of the tip. In the case of the 
power dependence 

1(p) = WW 

we get 

2ttR C 2nRd . . 

- 7 (d) 



n — 1 d n 1 n — 1 

In a more general case one must use numerical integration. 

In the local optic approximation, where the dielectric function is assumed 
to depend only on the frequency uj, the reflection factors R p and R s for flat 
surfaces, covered by an adsorbate layer, are given by [26J: 

p- s/e- A-Kin a q[sa\\/t- qa ± ] 

ttp — ; ; r, I ZD I 

" - 1 - s/e — Amn a q[saye + qa±\ 



P 



where 



R = P~ s -47rm a (^/c) 2 a|| 
s p + s + 4irin a (uj / 'c) 2 ct|| ' 



(28) 




and where aii and a± are the polarizabilities of adsorbates in a direction 
parallel and normal to the surface, respectively. Here e = e(uj) is the bulk 
dielectric function and n a is the concentration of adsorbates. For clean sur- 
faces n a = 0, and in this case formulas (??) reduce to the well-known Fresnel 
formula. 

At d < l,vp/u> and kp ~ 1, where I is the electron mean free path, and 
where Vp and kp are the Fermi velocity and Fermi wave number, respectively, 
the system will be characterized by non-local dielectric function e(q, to). In 
this paper we use the non-local optic dielectric approach, proposed some 
years ago for the investigations of the optical properties of a semi-infinite 
electron gas |2T)] . 

Accordingly to [20 , the reflection factor for p - polarized electromagnetic 
field, incident on the flat surface, is determined by (201 

* = wk- (29) 
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where the surface impedance Z p is given by 



= 2i r°° dq, ( q 2 {u/c) 2 g 2 z » 

^" vr 7o Q 2 \*i(u,Q) + (cj/c) 2 ei(cj,Q) -Q 2 1 l ° ' 

where t\ is the finite- life- time generalization of the longitudinal Lindhard 
dielectric function which accordingly to |2Sj can be written as: 

A ' ^ 1 + + i/r, g) - l)/(e°(0, Q) - 1) ' W 



3^2 

6 Z V,g) = i + ^|^ (32) 



1 , 1 / \2i i z ~ u + 1 , n / \2n ^ + 

^ = o + F" 1 " z " M ln 7" +!-(* + w ln — 7 

2 8,2 V 2 — -u— 1 z + w — 1 

■)2 _ J2 i „2 



(33) 

where Q 2 = q 2 + q 2 ,z = Q/2k F ,u = u/(Qvp), u p is the plasma frequency, 
r is the Drude relaxation time, and where vf and fcp are the Fermi velocity 
and wave vector, respectively. For s— polarization the reflection factor is 
determined by 

*■ = < 34 > 

1 + Z s p 

where 

y 2z f°° dq z 

s " vr Jo (u/c) 2 e t (uj,Q)-Q 2 { ) 

e*(w,g) = 1- , ^ . r /t, (36) 
a; (a; + 17) 

+[1 - (2 + „ 0Tln i±^±i) ,37, 

with -u' = (u; + ir _1 ) / (Qvp)- We show in Sec. 4 that the maximum of the elec- 
tromagnetic friction is reached for small electron densities, where the electron 
gas becomes non-degenerate (the electro gas is degenerate for ksT « Ef 
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and non-degenerate for fcgT > ep, where Bp is the Fermi energy). For non- 
degenerate electron gas we use the following classical expressions for dielectric 
functions [21] 

et(w,Q) 

where the function F(x) is defined by the integral 

F(x) = -4= / dz— (40) 

and Vt = JkBT/m, where m is the electron mass. 



1 + 



to. 



2 r 



Qv T , 



1 + F 



\V2Qv T , 

F ( ^ 



LO[LO 



vt 



(38) 
(39) 



4 The case of the good conductors 

By a well-conducting metal we mean one whose dielectric function e = 1 — 
4nia/uj (a is the conductivity) has an absolute value much larger than unity. 
For good conductors at thermal frequencies R P i « 1 and R pr ~ 1. Thus an 
enhancement in friction is possible only for very small q « 1/d . 

It is convenient to write the friction coefficient for the two flat surfaces in 
the form 

-" = »f *" (-£)<'' +I -> (41) 



Taking into account that qdq = kdk, from Eq. (J22J) for normal relative 
motion of clean surfaces within local optic approximation we get contribution 
to friction from evanescent p-and s-polarized electromagnetic waves 

POO rjh' r 

I17 n = / A; 5 [Re(s/e)] 2 [[(A; 2 + |s/e| 2 )cosh/W + 2A;[Im(s/e)]sinhA;d] 2 

j0 7T L 

+ (k 2 - \s/e\ 2 ) 2 } 5 (42) 

J |((s/e) 2 — k 2 )smhkd + 2i/c(s/e)cosh/c<i| 

ITs n = / -7 k 5 [Res} 2 [{k 2 + |s| 2 )cosh/W + 2fclmssinh£;d] 2 

JO IT L 
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+ (k 2 - \s\ 2 ) 2 j (43) 

J |(s 2 — k 2 )sinhkd + 2ikscoshkd\ 

As k — > 0, there is no singularity in I±^ n , and for I±p n , for (c/^)! 6 ! -3 ^ 2 < 
d < (c/ o^) lei 1 / 2 , than taking into account that in this limit sinhkd ps A;rf and 
cosh/crf ps 1, we get 

/•°° 



^(I+^^'-ttwf)' (44) 

where the surface impedance £ = e -1 / 2 = £' — For d < (c/u;)|e| -1 / 2 , 7 S 
becomes slowly dependent on d: 

/± s » -^(^/c) 4 | e | 2 (1.22 - ln(2d|e| 1 /Vc) (45) 

For d > (c/o;) |e| -1 / 2 we get 



/± s « (c/c) 2 C' 2 rf- 6 (46) 



For the propagating electromagnetic waves, taking into account that qdq = 
—pdp, we get 

IT* = (ulcfC 2 ^ % 5 l + cos 2 (pd) 

±p 1 / M 7o 7r 2 ^ |psinj9d + 22(cc;/c)Ccosprf| 4 ' 1 J 

rrad _ /, , / rW >2 [ U/C dp 5 1 + COS 2 ( P rf) 

±s l /Cj Wo 7T 2:P |(cu/c)sinpd + 2zpCcosp(i| 4 1 J 

For d < (c/u;)|e| -1 / 2 the contribution to friction from propagating wave is 
negligibly small in the comparison with the contribution from the evanescent 
waves. For d > (c/u;) |e| -1 / 2 as p — > 0, the integral I™ s d has no singularity, 
and we get for I* 



rad 



V«^{h^™ + TT$?r)- (49) 
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In addition, I rad has singularities at the other zeroes of sm.pl, i.e., near the 
values p n = nn/d < w/c (n is an integer). In the vicinity of p n , putting 
P = Pn + V ' i we have smpd ~ (— 1)" and cospd pa (— l) n , 

rrad _ 9 / I \ 2 r l2 f <¥_ 1 

~ f ![ _ ar ctan C'VC' Cto^ I (50) 

8tt 2 ^C v 2 1 + (C70 / 

The number m of such contribution is obviously equal to the integer part 

of the quantity y = uod/ixc (m = [y]), so that all p n (with the exception of 

p — 0) make a summary contribution 

! — arctan C'VC ~~ ~ ' 



x 



8cud%'\2 s/s 1 + (C70V 

'{m+lf (m + l) 4 (m+ l) 3 m 1 



5 2 3 30 30 

In the integral I^f, the contribution from the vicinity of the point p n is 

J /° dp' p\ 



(51) 



2(c/c)V 2 / ^ 

JO 7T 



2 \(u/c)p'd + 2ip n (\ 4 
-, I - - arctan C'VC 



8c#c'V 2 b/b i + (C70V 

and consequently 

= 8^? fi" arctanC7C ' - T+Sp) ? " 2 = 

" -arctan C7C'- ^—77^777^ ] m(m + l) (2m + 1) (52) 



48c#C' I 2 s/s 1 + (C'VC) 2 / 

At m » 1, when nuj/cd, the s— and p— wave contribu- 

tion are approximately equal, and for the total contribution from propagating 
electromagnetic waves in this limit we get 

1 1 LJ 4 

Trad rrad , rrad ^ 

J± ~ ^ + ^ s ~ 240vr3 c 4 C ' (53) 
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The above formulas in this Section were obtained from the general Eqs. 
(|21H24|) assuming absence of spatial dispersion of the dielectric function. But 
these formulas contain only surface impedance ( that describe the ratio of the 
tangential components of the electric and magnetic fields on the boundary of 
the body. Thus, the results in this section remain in force also in the presence 
of spatial dispersion, provided only that the surface impedance of the medium 
is small enough. We would have arrived at the same formulas also if we were 
to assume from the very beginning that the Leontovich boundary condition 
E = (H x n is satisfied on the surface of the metal. 

At not too low temperatures, the impedances of metals are given by the 
formula for normal skin effect 

C = C" = (uj/Sna) 1 / 2 (54) 

where a is the conductivity. In the local optic approximation it is assuming 
that there is no dependence of a on q. In the Wien region of frequencies it is 
also good approximation to neglect by the frequency dependence of a. In this 
approximation using flUJ) for \ w (k B T/4Tcha) 3/2 < d < \ w (4nh<j/(k B T) 1/2 
(Xw = ch/(ksT) we get 

^ = fi f ("I) - 0.13^ g£) ,55) 

For the comparison the p-wave contribution for parallel relative motion for 

is given by ([HUE!) 

It is interesting to note that for normal relative motion in contrast to the 
parallel relative motion practically for all d > the main contribution to 
friction comes from retardation effects because Eqs. in contrast to Eq. 
(j5Tlj) contains the light velocity. 

From Eq. (|4"5|) we get s-wave contribution to friction for d < X c 

7 evan ^ 10 -^( 3 _ 5 ln(2d/A c )) (57) 

For parallel relative motion the s-wave contribution is in two times smaller. 



16 



For d > A c , taking into account that Eq. (J46|) is valid only for uo > 

c 2 /An ad 2 , we get 

evan ^ B To 

ls d*C* ^ 

From Eq. ()53j) we get distance independent contribution to friction from 
propagating electromagnetic waves for d > X w 

Y p a J p ^1.9-10- 2 -^- (59) 

Figures 1-2 show the calculated contribution to the friction coefficient 
7 from evanescent electromagnetic waves for two semi-infinite solids, with 
parameters chosen to correspond to copper (r _1 = 2.5 • 10 13 s _1 , u p = 
1.6 • 10 16 s _1 ) at T = 273 K, for parallel (Fig.l) and normal (Fig. 2) relative 
motion. Results are shown separately for both the s- and p- wave contri- 
bution. The dashed line show the result when the local (long-wavelength) 
dielectric function e{uj) = t\ = e t is used, where 

u 2 

e(u) = 1 - — (60) 

In this case the integration in (|30p and (|35|) can be performed analytically 
and we get Fresnel formulas. Fig. 1 shows that the non-local optic effects 
become important for parallel relative motion for the p— wave contribution 
for sufficiently small separations (d < 1000 A). However, for the s— wave 
contribution for both parallel and normal motion the non-local optic effects 
are negligibly small for practically all separations. For normal relative motion 
for p-wave contribution the non-local optic effects are less important than 
for the parallel relative motion. In the present calculations we use the non- 
local dielectric approach which take into account the non-local optic effects 
on the length scale of the skin-depth. There are also the short range non- 
local optic effects coming from the non-local nature of the screening response 
near the surface. This gives the surface non-local contribution which we 
investigated in our previous publications [TTJ H2] . Comparing our previous 
calculations with the present one, we find that for d > 10A the volume 
contribution from the non-local effects is of the same importance as the 
surface contribution. For d > 10 A the main contribution to the friction 
coefficient 7 comes from s-polarized waves. In particularly, at d — 100 A the 
s-wave contribution 7 S w I0~ 5 kgs~ 1 m~ 2 , so that with the surface area A w 
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10 _14 m~ 2 ( typical for probe scanning microscopy), the friction coefficient is 
T m 7 S A ~ 10~ 19 A;<7S _1 . The s— wave contribution is characterized by weak 
distance dependence for d < 100A, and 7 ~ d- 2 for d > 100A. For good 
conductors like copper, even for very short distances, the main contribution 
to the friction coefficient comes from the s— polarized electromagnetic waves. 
This difference between p— and s— polarized waves results from screening 
effects: Good conductors are good reflectors for p— polarized field, which 
implies that they are ineffective in the emission and absorption of evanescent 
p— polarized waves. However these screening effects are less important for 
s— polarized waves. 

As pointed out in piJJ O E3 123, the p -wave contribution increase and 
the s— wave contribution decrease when the free electron density decrease. 
Within the local optic approximation the force of friction diverges when one 
go to the limit of zero conductivity. This situation is different from the radia- 
tive heat transfer, where, even in the case of local optics, a maximum in the 
heat transfer occurs for conductivities corresponding to semi-metals. Fig. 3 
shows the dependence of coefficient of friction on the electron density. When 
the electron density decreases there is transition from a degenerate electron 
gas to a non-degenerate electron gas at the density np ~ (fceTm) 3 / 2 /7r 2 h 3 . 
At T = 273 K the transition density np ~ 10 25 m~ 3 . For n > np we use 
the (non-local) dielectric function appropriate for a degenerate electron gas, 
while for n < np we use an expression corresponding to a non- degener- 
ate electron gas. In the calculations we used the electron mean free path 
/ m 600 A. At d — 100A the maximum value 7 maa . ~ 10~ 4 kg ■ s _1 is obtained 
for n max ~ 10 22 m~ 3 , corresponding to the DC conductivity o ~ \(Vt ■ m)~ l . 

5 Photon tunneling enhancement of the van 
der Waals friction 

We rewrite the denominator in Eq. ()22)1 in the form 

|1 _ e -^Rf = [(1 _ e - kd R r f + e~ 2kd R 2 ] 2 

x [(1 + e- kd R r f + e~ 2kd R 2 } 2 (61) 

where R r and Ri are real and imaginary part of R, respectively (R = R r + 
iRi). Let us suppose that \R r \ » Ri. In this case resonant conditions are 
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determined by equation 

R r (u±(q) = ±e qd (62) 
Close to resonance we can write 

[(1 ± e~ qd R r ) 2 ± e~ 2qd R 2 ] 

« e- 2qd R' r 2 (u ± )[(u - uj ± ) 2 + {R^/Ki^)) 2 } (63) 

where 

n/ . , dR'Ju) 
R r [uJ±) = 



du 

U=UJ± 

which leads to the following contribution to the friction coefficient: 

fi 2 rq c p 2kd 

^ 16k B TJo {\R' r (u±)\R i (u±)smh 2 hu±/2k B T] V ; 

The parameter q c in this expression defines the value of q where the two 
poles approximation is valid. To proceed further let us make the following 
simplifications. Close to the poles we can use approximation 

R = (65) 

UJ — UJq — IT] 

where a is a constant. Then from resonant condition (|6*2*|) we get 

uj± = lo ± ae~ qd 

For the two poles approximation to be valid the difference Aou = \u + — oj_\ 
must be greater than the width of the resonance 77. From this condition we 
get q c < \n(2a/rj)/d. For the short distances the parameter q c defines the 
value of q where the solution of Eq. (|62j) ceases to exit. 
For ujq > a and q c d > 1 from Eq. (jfHj) we get 

3 h 2 a 2 1 
7±± 128d i k B Tr ] smh 2 (htu /2k B T) { ' 

For the parallel relative motion using the same approximation we get 

7,1 128nk B Tsmh 2 (huj /2k B T) 1 ' 
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Interesting, the explicit d dependence has dropped out of Eq. JH7J). How- 
ever, 7|| is still d dependent, through the d dependence of q c . For the small 
distances one can expect that q c is determined by the dielectric properties of 
the material and does not depend on d. In this case the friction will be also 
distance independent. Probably, the weak distance dependence observed in 
[U] can be explained by the resonant photon tunneling. 

Resonant photon tunneling enhancement of the van der Waals friction 
is possible for two semiconductor surfaces which can support low-frequency 
surface plasmon modes. As an example we consider two clean surfaces of 
silicon carbide (SiC). The optical properties of this material can be described 
using an oscillator model [2H] 



rU">: r x (l+ / L f T . v ) (08) 
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s 



with eoo = 6.7, u L = 1.8 ■ lO 14 *" 1 , u T = 1.49 • 10 14 s-\ and T = 8.9 • 10 
The frequency of surface plasmons is determined by condition e r (u p ) = — 1 
and from (JSJ we get uj p = 1.78-10 14 s _1 . In Fig. 2 we plot the friction coefficient 
j(d): note that the friction between the two semiconductor surfaces is several 
order of magnitude larger than between two clean good conductor surfaces. 

Another enhancement mechanism is connected with resonant photon tun- 
neling between adsorbate vibrational modes localized on different surfaces. 
As an example, let us consider ions with charge e* adsorbed on metal sur- 
faces. The polarizability for ion vibration normal to the surface is given 
by 

e* 2 

a± = M^i-cu 2 -^)' (69) 

where u± is the frequency of the normal adsorbate vibration, and r] ± is the 
damping constant. In Eq. (J26|) the contribution from parallel vibrations 
is reduced by the small factor 1/e. However, the contribution of parallel 
vibrations to the van der Waals friction can nevertheless be important due to 
the indirect interaction of parallel adsorbate vibration with the electric field, 
via the metal conduction electron [27|. Thus, the small parallel component 
of the electric field will induce a strong electric current in the metal. The 
drag force between the electron flow and adsorbates can induce adsorbate 
vibrations parallel to the surface. This gives the polarizability: 

oq| = 7 2 \ ( 70 ) 

n e (oof, — u z — tiorju) 
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where n is the conduction electron concentration. As an illustration, in 
Fig. 3 we show coefficient of friction for the two Cu(OOl) surfaces covered 
by a low concentration of potassium atoms ( n a = 10 18 m~ 2 ) . In the q— 
integral in Eqs. (j22l24j) we used the cut off q c ~ it /a (where a ~ lnm is 
the inter-adsorbate distance) because our microscopic approach is applicable 
only when the wave length of the electromagnetic field is larger than double 
average distance between the adsorbates. In comparison, the friction between 
two clean surface at separation d = lnm is seven order of magnitude smaller. 
At d = lnm the friction coefficient T for an atomic force microscope tip with 
R ~ 1/im is ~ 10~ l2 kgs~ l (7 ~ 10 3 kgs~ l m~ 2 , see Fig.2); this is of the same 
order of magnitude as the observed friction jHj. 



6 Summary and conclusion 

We have calculated the van der Waals friction between two flat surfaces 
for normal relative motion and found a drastic difference in the comparison 
with parallel relative motion. This difference is connected with resonance 
produced by the multiple scattering of the electromagnetic waves from the 
opposite surfaces. In the case of sharp resonance it gives much larger con- 
tribution to friction in the case of normal relative motion than for parallel 
relative motion. 

We have studied in the detail the friction between two good conductors 
and have found that for normal relative motion even for very small distances 
the main contribution to friction comes from the retardation effects. The 
non-local optic effects are very important for p-wave contribution to friction 
for parallel relative and much less important for normal relative motion. 
For s-wave contribution the non-local optic effects are unimportant for both 
direction of relative motion. 

In the case of van der Waals friction we have found that for distances 
between bodies ~ 100A, for good conductors with a high free electrons 
concentration, the main contribution to friction is associated with the s- 
polarized electromagnetic waves. For d < 100 A this mechanism gives a 
friction coefficient per unit area 7 ~ 10 _5 kgs _1 m~ 2 , nearly independent of 
the distance d, while for d > 100A the friction coefficient 7 depends on dis- 
tance as d~ 2 . For an atomic force microscope tip with the near substrate 
area A ~ 10 _14 m 2 , we got the friction coefficient T ~ 7 A ~ 10 _19 kgs _1 for 
d < IOOA. When the concentration of electrons decreases, the s— contribution 
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to the friction decreases while the p— contribution increases. At d = 100A 
and with the electron lifetime r = 4 • 10 _14 s, the p— contribution reaches 
maximum 7 max ~ 10~ 4 kgs _1 m~ 2 at the electron concentration n ~ 10 22 m~ 3 , 
which corresponds to the conductivity a ~ 1 (fim) -1 . 

We have shown that the van der Waals friction can be enhanced by sev- 
eral orders of magnitude in the case of resonant photon tunneling between 
low-frequency surface plasmon modes and adsorbate vibrational modes. In 
the case of friction for two Cu(100) surfaces covered by a low concentration 
of potassium atoms at d — 10A we have found friction of the same order 
of magnitude as it was observed in experiment jHj. However, the distance 
dependence in this case is more stronger than it was observed in [U] . Further 
experiments with simple and well defined composition of the tip and sample 
must be performed to elucidate different energy dissipation mechanisms in 
the non-contact friction. 

The obtained results should have broad application in non-contact friction 
microscopy, and in design of new tools for studying adsorbate vibrational 
dynamics and optical properties of the surface plasmons. 
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After substituting J3J) and (@J) into formula (|19|) to linear order in vibrational 
coordinate u and frequency u> we get 



47T JO 



du 



+ < \Voz\ >) 



UJ 



+ (p -/>*)(< w 0z v* 0z > +c.c)]+ - ) p (p + p*)(< \woy\ 2 > + < \v 0y \ z >) 



i f p r 

+{P-P*){< WoyV^ > +c.c)J + [-^[(p + p*)(< W 1z Wq z > + < V 1z Vq z > +C.C.) 
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(P-P*)(< W 1zVq z > + < V 1z Wq z > +C.C.) 



-p 



{P+P*)(< WlyW* Qy > 

—iujot 



lu(lu + LO ) 

+ < V 1z Vq z > +C.C.) +{p~ p*)(< W ly VQy > + < V ly wl y > +C.C 

(71) 

From Eqs. (J2*U|) and ()71|) it follow that the friction coefficient is determined 
by the formula 



1± 



d 2 q t d (p 



(^7;(^[(P+P*)(< w Xz wl z > + < v lz v* 0z > -c.c.) 

(P+P*)(< WlyWQy > 



;P 



Anuoi Jo "J (2n) 2 ^dcuo ^ q 

(p-p*)(< %i > + < > ~ c - c v 

+ < v lz v* 0z > -c.c.) + (p-p*)(< w ly v* 0y > + < v ly w* 0y > -c.c.)])) (72) 
Using Eqs. (|11II15I17I18|) we get 
1 g 
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(p+(p+p*)(< w lz w* 0z > + < v 1z Vq z >) - c.c) wo=o = 2ra (n(w) + 



(1 - 


RipR2p 


2 ) 2 + l(l- 


-Rip 


2 )i? 2p e^ + (1 - 


-R2p 


2 )ii!^e- ipd | 2 




Ap 


4 



x 77~P 



g 2 <9co»c 



w =0 



X 



9 p 2 

aJ|A p l 4 



(Imi2 lp + e- 2|p|c! | J Ri p | 2 Im J R 2p )(Im J R 2p + e - 2Md \R lp \ 2 ImR 2p ) 



+ e- 2 ^ d lm(R lp R 2p ) 2 \e- 2 ^ d (74) 

Other similar expressions for the s-wave contribution can be obtained from 
Eqs. (|73^) and (J71J) by replacement of the reflection factors R p for p-polarized 
wave by the reflection factors R s for s-polarized wave. After substituting 
Eqs. (J73J) and (fT^) . and similar expression for s-polarized waves in Eq. (J72|) 
we get the friction coefficient for normal relative motion which is determined 
by formulas (ET H22jl . 



FIGURE CAPTIONS 
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Fig. 1. The friction coefficient for two flat surfaces in parallel relative 
motion as a function of separation d at T = 273 K with parameter cho- 
sen to correspond to copper (r~ l = 2.5 ■ 10 13 s _1 , uj p = 1.6 ■ 10 16 s _1 ). The 
contributions from the s— and p— polarized electromagnetic field are shown 
separately. The full curves represent the results obtained within the non- 
local optic dielectric formalism, and the dashed curves represent the result 
obtained within local optic approximation. (The log-function is with basis 
10) 

Fig. 2. The friction coefficient for two flat surfaces in normal relative 
motion as a function of separation d at T = 273 K with parameter cho- 
sen to correspond to copper (r _1 = 2.5 • 10 13 s _1 , u p = 1.6 • 10 16 s _1 ). The 
contributions from the s— and p— polarized electromagnetic field are shown 
separately. The full curves represent the results obtained within the non- 
local optic dielectric formalism, and the dashed curves represent the result 
obtained within local optic approximation. (The log-function is with basis 
10) 

Fig. 3. The friction coefficient for two flat surface as a function of the 
free electron density n at T = 273 K. The full curve was obtained by in- 
terpolation between the result (dashed lines) obtained within the non-local 
optic dielectric approach, with dielectric functions corresponding to a degen- 
erate electron gas for n > rip ~ 10 25 m -3 , and to a non-degenerate electron 
gas for n < rip. The calculation were performed with the damping constant 
t- 1 = 2.5 • lO 13 ^ 1 , separation d = 100A and n = 8.6 • 10 28 m" 3 . (The 
log-function is with basis 10) 

Fig. 4. The friction coefficient for two clean semiconductor surfaces in (a) 
normal and (b) parallel relative motion, as a function of the separation d. 
T = 300K and with parameters chosen to correspond to a surfaces of silicon 
carbide (SiC) (see text for explanation) (The log-function is with basis 10) 

Fig. 5. The friction coefficient for two surface covered by adsorbates in 
(a) normal and (b) parallel relative motion, as a function of the separation 
d. T = 273K and with parameters chosen to correspond to K/Cu(001) [29 
(u ± = 1.9- 10 13 s-\u;|| = 4.5-10 12 s- 1 ,r/|| = 2.8-10 10 s~ 1 ,r ]± = l.Q-10 12 s~\ e* = 
0.88e) (The log-function is with basis 10) 
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